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We present studies of the potential energy landscape of selected binary Lennard- Jones thirteen 
atom clusters. The effect of adding selected impurity atoms to a homogeneous cluster is explored. 
We analyze the energy landscapes of the studied systems using disconnect ivity graphs. The required 
inherent structures and transition states for the construction of disconnectivity graphs are found by 
combination of conjugate gradient and eigenvector- following methods. We show that it is possible 
to controllably induce new structures as well as reorder and stabilize existing structures that are 
characteristic of higher-lying minima. Moreover, it is shown that the selected structures can have 
experimentally relevant lifetimes. 

PACS numbers: 



I. INTRODUCTION 

The science of chemistry is characterized by an inter- 
play of reductionist and constructionist themes. On the 
one hand, one seeks to reduce complex systems to more 
understandable and more controllable components. On 
the other, one strives to utilize these components to con- 
struct complex assemblies to meet specific chemical, bi- 
ological, and/or materials goals. The core elements of 
such efforts have, in the past, been largely atomic and/or 
molecular in nature. 

Viewed from the above perspective, a central compo- 
nent of much ongoing research has been the develop- 
ment of a variety of "reductionist" classical and quantum- 
mechanical tools. Common goals in the application of 
these tools have been locating the minimum energy con- 
figuration for the associated potential energy surface^ 
and sampling the relevant, low- lying local minima.^i^i^ 

Increasingly, both chemistry and materials science are 
entering an era in which the fundamental components of 
the constructionist phase of the problem often are them- 
selves complex, preassembled objects. Sun, et al.^ioi ex- 
ample, have shown that novel magnetic materials can be 
prepared through the self-assembly of colloidal clusters. 
In related developments Lehmann and Scoles^ as well 
as Miller, et ai^'^'^^ have demonstrated the effectiveness 
of superfluid solvents in preparing unusual, metastable 
species. Sustained progress in this emerging field of clus- 
ter assembled materials ultimately rests on the ability 
to characterize and control a broad range of increasingly 
complex nanoscale objects. 

In a previous paper^^, we have examined a number 
of theoretical issues of general concern with respect to 
predicting/characterizing/controlling the structure and 



dynamics of cluster-based precursors. That work, in 
essence, seeks to invert the logic of the minimization 
problem. That is, instead of searching for the mini- 
mum energy structures of specified energy landscapes, we 
strive instead to reshape those landscapes and thereby to 
exercise control over selected physical systems. In par- 
ticular, we seek to stabilize and/or kinetically trap con- 
formers of the parent homogeneous system that are oth- 
erwise either un- or metastable. Our hope is that such 
new structures might have interesting physical proper- 
ties (electronic, magnetic, optical, thermal, etc.) and 
can be used as precursors in subsequent assembly pro- 
cedures. Our previous work^^ has demonstrated that for 
small Lennard- Jones systems we could, using selected im- 
purities, both alter the energy ordering of the stable core 
atom isomers and induce wholly new conformers not seen 
in the original homogeneous species. We have also seen 
similar reordering effects in our studies of molecular ni- 
trogen adsorbates on nickel clusters. 

In the present study, we wish to extend the results of 
our previous investigations. First, we wish to consider a 
number of larger clusters to demonstrate the general ap- 
plicability of our efforts. Second, we wish to demonstrate 
the resulting species can be made sufficiently robust that 
they are of practical interest. To attain the second goal 
we enhance our previous publication^^ by determining 
RRKM isomerization rate constants and lifetimes for the 
generated isomers. 

An important byproduct of the current work is the in- 
formation we garner concerning the structure of the po- 
tential energy surfaces of the mixed cluster systems. In 
the current work we find that many of the disconnectiv- 
ity graphs have double- funnel structures. The thermo- 
dynamic properties of other cluster systems with double- 
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funnel structureaSii^ can be rich exhibiting sohd-hke to 
sohd-hke phase change phenomena. Motivated by these 
past studies and the disconnectivity graphs investigated 
in the current work, in a companion paper^^ we investi- 
gate the energy and heat capacity of some of the systems 
studied here. 

The remainder of the paper is organized as fohows. 
An outhne of the computational details of this work is 
presented in Sec^ We discuss the methods utilized to 
find the inherent structures and transition states on the 
potential energy surface as well as the method to estimate 
the inherent structures lifetimes. In Sec. IIIII we present 
the results that demonstrate "proof of principle" with 
respect to the goals of the present studies. In Sec. II VI we 
summarize our results and speculate about likely future 
research directions. 



II. COMPUTATIONAL DETAILS 

In this section, we describe the computational de- 
tails of the studies involving binary clusters of the form 
Xi3_nYn. Our overall interest is to explore the extent to 
which we can utilize the "adatoms" (i.e. the Y-system) to 
induce, reorder and stabilize selected inherent structures 
in the "core" X-system. While one can easily imagine 
applications involving more and more complex compo- 
nents, we have found that these relatively simple, two- 
component clusters are a convenient starting point for 
an initial study of the issues we raise. 

The total potential energy, Vtot^ of a cluster consisting 
of N particles is modeled as a pairwise sum of Lennard- 
Jones interactions 

N 

Vtot = '^Vij{rij), (1) 

i<j 

where the pair interaction as a function of the distance 
between particles i and j, r^j, is given by 

%-(r'.,) = 46,, [(^)^^-(^n. (2) 

In Eq. ((2j the constants e^j and cr^j are the energy and 
length-scale parameters for the interaction of particles i 
and j. 

For a binary system, both the "like" (X-X, Y-Y) as 
well as the "unlike" (X-Y) interactions have to be spec- 
ified. With an eye toward studying trends in the results 
as opposed to results for particular physical systems, it is 
convenient to reduce the number of free parameters. To 
do so, we shall assume in the present study that the "un- 
like" Lennard- Jones values are obtained from the "like" 
Lennard- Jones parameters via usual combination rule&i^ 

crxY = -^{o'xx + cryy), (3) 



^XY = ^/^xxJyY' (4) 

Furthermore, we note that with the mixed Lennard- Jones 
parameters specified as in Eqs.(|3|) and Q), the inherent 
structure topography of the "reduced" potential energy 
surface of the binary system (i.e. Vtotjexx) is a function 
of only two parameters, (cr, e), the ratios of the corre- 
sponding adatom/core length and energy parameters 

(7-=gyyIgxx, (5) 



e = eyy/exx. (6) 

If necessary for a discussion of a specific physical sys- 
tem, the absolute bond lengths, energies, activation en- 
ergies, etc. can be obtained from the corresponding "re- 
duced" results by a simple rescaling with the appropriate 
core-system Lennard- Jones parameters. 

A. Stationary points and the disconnectivity graph 

The computational task in our study is thus one of ex- 
ploring and characterizing the (reduced) potential energy 
surface of our binary cluster systems as a function of the 
number of (core, adatom) particles, (n, m), and for given 
(cr, e) ratios. In typical applications the lowest N/^ inher- 
ent structures and the associated disconnectivity graphs 
are determined. For the applications reported here, N/^ 
is generally of the order of several thousands or more. 
The inherent structures are found either via conjugate 
gradient methods^ starting from randomly chosen ini- 
tial configurations, or by more systematic surface explo- 
ration methods. ^^'^^ In all cases, the inherent structures 
that are located are confirmed to be stable minima via a 
standard Hessian analysis. To reduce the chance we miss 
particular local or global minima, we monitor the num- 
ber of times individual inherent structures are found and 
demand that each of the N/^ inherent structures be lo- 
cated a minimum number of times (at least 10) before we 
terminate our search. Once we are satisfied we have lo- 
cated the relevant inherent structures, transitions states 
linking these stable minima are obtained using the eigen- 
vector following methods outlined by Cerjan and Miller^^f 
and further developed by Simons et alS^^^, Jordan et 
alw^ and Wales^^. Finally, with the given inherent struc- 
tures and transition states, we perform a disconnectivity 
analysis 

B. Rate constants and lifetimes of the inherent 
structures 

From the known inherent structures and the transition 
states that connect /separate them, we estimate rates for 
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transitions between neighboring inherent structures. The 
rate constants allow us to calculate the average amount of 
time the system will spend in a given inherent structure, 
i.e. the lifetime of an inherent structure. 

There are variety of methods available to estimate 
the rate constants (see R,efl26l and references therein). 
We utilize the harmonic approximation to the Rice- 
Ramsperger-Kassel-Marcus (RRKM) method. It has 
been found that this method gives good estimates of rates 
for isomer izat ion of clusters^^. The rate constant, /c^j, for 
transition leading from inherent structure j to inherent 
structure i is given as a sum over all transition states 
connecting inherent structures i and 

ki,=J2kf. (7) 

a 

is given by 

la i-rSAT-e IS 
Ul=l 

where A(j)j is 

A^j = Ef - Ef , (9) 

Ej^^'" is the energy of the transition state, E?^ is the 
energy of the inherent structure, and uij are the cor- 
responding normal mode frequencies. N is the number 
of particles in the cluster while hj and are the order 
of the point group of inherent structure j and transi- 
tion state a, respectively. Since we are interested in the 
qualitative estimates of the rate constants (lifetimes) we 
neglect the order of the point group of inherent struc- 
ture and transition state from our calculations. We es- 
timate that the error arising from this approximation is 
a factor between 60 and 1 based on the reasoning that 
follows: The global minimum of the systems we study 
is icosahedral in nature which means that the order of 
its point group is no larger than hj =60. The transition 
state that connects the global minimum with a higher 
lying inherent structure has a lower symmetry than the 
global minimum and therefore smaller order of the point 
group than the one associated with the point group of the 
global minimum. As is made evident below, the factor 
of 60 is unimportant to the determination of the order 
of magnitude estimates of isomer lifetimes that are of 
interest to us in the current context. We would like to 
point out that the numerical value of exx and axx used 
to calculate the rate constants are 119.8 K and 3.405 A, 
respectively. 

III. NUMERICAL RESULTS 

In the present Section, we illustrate the general themes 
we introduced in Section^. We demonstrate that we can 
accomplish three basic objectives. Specifically, we show 



that by adding selected "impurity" atoms to bare "core" 
systems, we can: 

1. induce new "core structures" 

2. reorder the energies of existing core inherent struc- 
tures, and 

3. stabilize selected inherent structures by controlling 
the activation energies that determine their isomerization 
kinetics. 

For purposes of illustration, we examine numerical re- 
sults for three, thirteen atom Lennard- Jones systems in- 
volving ten, eleven and twelve core atoms, systems well- 
known from previous studies^^ to have 64, 170 and 515 
energetically distinct inherent structures, respectively. 
These systems have been chosen because they build upon 
simple ten, eleven and twelve-atom cores and because 
total systems have thirteen atoms, a magic number for 
icosahedral growth in homogeneous systems. 



A. X12Y1 

We first consider binary clusters of the type X12Y1. 
The selected inherent structures and their associated en- 
ergies for X12 core system are illustrated in Fig. ^ The 
inherent structure labeled by (a) is the global minimum 
while all others are the higher lying inherent structures. 
Here one impurity atom Y is added to the parent, twelve- 
atom X core. Using the techniques of Section^ we de- 
termine the lowest several inherent structures for a range 
of (cr,e) [c.f. Eq.(|S|) and Eq.©]. As can be seen from 
Fig. [3 the total potential energy [Eq.(P)] of the lowest 
inherent structure for the X12Y1 system shows no ap- 
preciable structure as a function of the (cr,e) parameters. 

On the other hand, it can be seen in Fig.|31that the core 
potential energy, defined as the potential energy of inter- 
action for only the core X-atoms, of the minimum (total) 
energy cluster clearly breaks into extended regions, each 
corresponding to a well-defined core structure. We would 
like to point out that each region in Fig. |31 contains the 
same "kind" of core structure but their core energies are 
slightly different. We have chosen a single "average" core 
energy value to represent all energies in the correspond- 
ing domain for plotting convenience. 

The distinct core structures, shown in Fig. [31 have been 
identified by examining their core energies (Ecore) and 
their principal moments of inertia. For each structure 
a triplet of values (Ecore, I2, I3) has been associated, 
where I2 and I3 are the moments of inertia about the 
principal axes 2 and 3, respectively. We have defined I2 
and I3 in the following way: I2=l2/Ii, 13=13/11 where I^, 
I2 and I3 are the principal moments of inertia obtained 
by diagonalizing the inertia tensor of the system. If the 
triplet of values has not been sufficient to identify a core 
structure then we have examined the structure visually. 

Selected cluster structures illustrating the core ar- 
rangements corresponding to various (cr,e) values are 
shown in Fig. ^ It can be seen from Figs. El and ^ that 
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(a) (b) 



(c) (d) 




(e) (f) 



FIG. 1: The selected stable inherent structure for X12 LJ 
cluster. Their energies (in units of the LJ well depth) are: (a) 
-37.968, (b) -36.347, (c) -36.209, (d) -35.461, (e) -34.447, (f) 
-33.598. 




FIG. 2: Etot(cr,e) (c.f.EqsC] Eland for the X12Y1 system. 
Note the relative lack of structure in the (cr,e) variation of the 
total cluster energy. 




FIG. 3: Ecore('7, e) for the X12Y1 system. Here the "core" en- 
ergy is defined as that portion of the potential energy arising 
from only the core-core atom interactions. Unlike the total 
energy, the (a, e) variation of the core cluster energy exhibits 
relatively well-defined regions. The labels of each of these re- 
gions in the figure correspond to the distinct core structures 
shown in Fig. 0] 



the X12Y1 cluster exhibits core X-atom structures that 
are higher lying minima in the parent X12 system. The 
core structures labeled by (4.2), (4.3), (4.4), and (4.6) 
can be recognized as structures labeled by (c), (d), (e), 
and (f) in Fig. ^ respectively. This illustrates that a 
suitable choice of the (a, e) parameters can controllably 
reorder the energies of existing core inherent structures. 
The structure labeled by (4.5) shows a newly induced 
core geometry not present as a stable minimum in the 
bare cluster. These two results demonstrate that we can 
accomplish objectives (1) and (2) stated earlier. 

Figure [Slrepresents the X12Y1 cluster at four points in 
Fig.Oldefined by the (a, e) coordinates (0.8,0.6), (0.8,1.0), 
(0.8,1.5), and (0.8,2.0). Here the pairs of coordinates cor- 
respond to (a), (b), (c), and (d) of Fig. [51 respectively. In 
other words, we keep value of cr=0.8 fixed, while increas- 
ing the value of e. The number of inherent structures 
available to the X12Y1 cluster varies from at least 4153 
in Fig. [51a to at least 3641 in Fig. [51d. Since we are 
primarily interested in energetically low-lying inherent 
structures we show only the lowest 200 inherent struc- 
tures. The global minimum of each system is labeled 
by the number 1 and contains as a recognizable compo- 
nent the core structure labeled by (4.6) (see Fig. EJ. In 
Fig. [51a the core structure (4.6) is connected to a group 
of 12 inherent structures by pathways whose energies do 
not exceed -39.4 (in units of exx)- The energies of these 
12 inherent structures are very close to each other and 
their values are from lowest -39.810 to highest -39.614. 
Their corresponding core structures are different from the 
core structure associated with the global minimum. This 
implies that barriers that connect the global minimum 
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(4.1) (4.2) 




(4.3) (4.4) 



(4.5) (4.6) 

FIG. 4: Plots of X12Y1 structures for selected (a, e) values. 
The decimal number for each figure denotes the corresponding 
(cr, e) domain in Fig. El 



with the inherent structures in the group are "relevant" 
barriers. The relevant barriers as those that connect in- 
herent structures that contain different core structures. 
By examining the disconnectivity graph on a finer energy 
scale we find that inherent structure 1 is connected, by 
the lowest isomerization barrier, to inherent structure 16. 
The numerical value of the lowest isomerization barrier 
is AEi^i6=2.623exx- Figures. [5Jb and[5Jc show that in- 
creasing the value of e increases the isomerization barriers 
that connect inherent structure 1 with a group of 12 and 
9 inherent structures, respectively. Similarly to system 
(a), the inherent structures associated with both groups 
contain core structures that are different from the core 
structure associated with inherent structure 1. For sys- 
tem (b) the lowest isomerization barrier that connects in- 
herent structure 4 with inherent structure 1 has a numeri- 
cal value AEi^4=2.960exx while for system (c) the lowest 
isomerization barrier connects inherent structure 1 with 
inherent structure 6 and its value is AEi^6=3.269exx- In 



Fig. [5Jd inherent structure 1 is connected to a group of 
9 inherent structures. The lowest isomerization barrier 
is AEi^3=3.509exx connecting inherent structure 1 with 
inherent structure 3. The double- funnel structure of the 
disconnectivity graphs are evident especially in Fig. [^la 
and Fig. [Sib where the minima that define the two sep- 
arated basins are so close in energy. The double-funnel 
structure of the potential energy surface is reflected in 
the classical heat capacity as discussed in the companion 
paper— Below we find similar double- funnel structures 
for XnY2. 

We estimate the rate constants (lifetimes) for four tem- 
peratures, 5, 10, 100, and 300 K, as a function of the 
height of the isomerization barriers (e). At the low tem- 
peratures (5 and 10 K) the studied systems become ex- 
tremely stable. By increasing the isomerization barrier 
between the global minimum and the first higher lying 
inherent structure, from AEi^i6=2.623exx (see Fig.Ela) 
to AEi^6=3.269exx (see Fig. \^d) the lifetime increases 
by nine and four orders of magnitude in the case of 5 
and 10 K, respectively. To be more specific, at 10 K, the 
lifetime increases from the order of seconds to the order 
of days. This is illustrated in Fig. IT^a. 

As illustrated in Fig. Eland Fig.Ela the barriers that 
determine the isomerization kinetics are sensitive to the 
((j,e) values and can thus be at least partially controlled. 
Therefore, we have created selected structures that have 
experimentally relevant lifetimes. These results are spe- 
cific demonstrations of goal (3) stated earlier. 



B. X11Y2 

As a second illustration, we consider mixed clusters 
of the type X11Y2. This system builds upon a parent, 
eleven- atom system known to exhibit a set of 170, en- 
ergetically distinct inherent structures'^. The selected 
core inherent structures and associated energies for the 
stable Xii inherent structures are presented in Fig.jSl In 
Fig. [3 a (<j, e) contour plot of the core-atom potential 
energies of the lowest total energy X11Y2 clusters, again 
reveals the presence of definite "core-phases". As il- 
lustrated in Fig. IHl some of these regions correspond to 
various core structures present in the parent Xn system 
while others correspond to new structures not seen in 
the original, single-component cluster. We can see from 
Figs. [3 to El that the impurity Y atoms provide us with 
significant control over the relative ordering of the core 
energies of the parent Xn system. Moreover, since we 
can manipulate the isomerization barriers in the X11Y2 
systems, we can at least partially stabilize clusters that 
exhibit selected core structures with respect to isomer- 
ization. This is illustrated in Fig. El 

Figures Ela - ^d represent the X11Y2 cluster at four 
points in Fig. [3 with XiiY2(cr, e) coordinates (0.8,0.5), 
(0.8,1.0), (0.8,1.5) and (0.8,2.0), respectively. The num- 
ber of inherent structures available to the X11Y2 cluster 
in all four cases is more than 6000. We show only the 





FIG. 5: Disconnectivity graph for X12Y1 (cr,e) values demonstrating that we can control barriers for the selected inherent 
structures. The energy scale is in units of exx- The (cr, e) values for panels (a-d) are (0.8,0.6), (0.8,1.0), (0.8,1.5), and (0.8,2.0), 
respectively. Only branches leading to the 200 lowest-energy minima are shown. 






(c) (d) 




(e) (f) 

FIG. 6: The selected stable inherent structure for Xn LJ 
cluster. Their energies (in units of the LJ well depth) are: (a) 
-31.766, (b) -31.9152, (c) -31.9146, (d) -31.775, (e) -31.615, 
(f) -31.036. 




0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 



a 

FIG. 7: Ecore{cr,e) for X11Y2. Format for the plot is the 
same as in Fig. El 



(8.4) (8.6) 




(8.8) (8.9) 

FIG. 8: Plots of selected X11Y2 structures for various (cr, e) 
values identified in Fig.|3 The number of the structures corre- 
spond to the regions labeled in Fig.|7| The core structures for 
the systems labeled by (8.6) and (8.8) are not stable energy 
structures of the bare Xn system. 



lowest 200 inherent structures. The global minimum of 
each system is labeled by the number 1 and contains 
as a recognizable component the core structure shown 
in Fig. ISlS. In Fig. [UJa the global minimum, the core 
structure labeled by (8.8) in Fig. [SI is linked to inher- 
ent structure 5. Its core structure is different from the 
one associated with the global minimum. The isomer- 
ization barrier between them is AEi^5 = 1.015exx- In 
Fig. IHlb inherent structure 1 is connected to inherent 
structures 3 and 4 whose energies are almost degener- 
ate. Both inherent structures 3 and 4 contain core struc- 
tures that are different from each other and from the 
one associated with inherent structure 1. The isomer- 
ization barriers between inherent structures 1 and 3 and 
1 and 4 are AEi,3=2.025exx and AEi,4=2.011exx, re- 
spectively. Figures [Hlc and Fig. I^ld show that further 
increasing of the value of e increases isomerization barri- 
ers that link inherent structure 1 with inherent structure 
2. Numerically, these barriers are AEi^2=2.265exx and 
AEi^2=2.463exx, respectively. Estimated lifetimes are 
shown in Fig. E|b as a function of e. It can be seen 




FIG. 9: Disconnectivity graph for X11Y2 (cr,e) values demonstrating that we can control barriers for the selected inherent 
structures. The energy scale is in units of exx- The (cr, e) values for panels (a-d) are (0.8,0.5), (0.8,1.0), (0.8,1.5) and (0.8,2.0), 
respectively. Only branches leading to the 200 lowest-energy minima are shown. 
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that the hfetimes increase by fourteen and seven orders 
of magnitude at 5 and lOK, respectively. 



C. X10Y3 

As a third ihustration, we consider mixed clusters of 
the type X10Y3. This system builds upon a parent, ten- 
atom, system known to exhibit a set of 64, energetically 
distinct inherent structures'^. The selected core inher- 
ent structures and associated energies for the stable Xio 
inherent structures are presented in Fig. ^21 A (a, e) 
contour plot of the core-atom potential energies of the 
lowest total energy Xio Y3 clusters is shown in Fig. ^2 




(c) (d) 




(e) (f) 

FIG. 10: The selected stable inherent structure for XioLJ 
cluster. Their energies (in units of the LJ well depth) are: (a) 
-28.422, (b) -27.556, (c) -27.214, (d) -26.772, (e) -26.698, (f) 
-26.695. 

As illustrated in Fig. some of the domains in Fig. 
correspond to the core structures present in the parent 
system Xio while others correspond to new structures 




FIG. 11: Ecore(cr, e) for X10Y3. Format for the plot is the 
same as m Figs. Eland □ 



not seen in original, single-component cluster. Specifi- 
cally, the domains 1, 2, 3, 4, and 5 correspond to the 
core structures labeled by (a), (c), (f), (e), and (d) in 
Fig. ^1 respectively. Therefore, the impurity atoms Y 
provide significant control over relative ordering of the 
core energies of the parent Xio system. FiglT^ illustrates 
that, by choosing the appropriate set of (a, e) values we 
can manipulate the isomerization barriers in the selected 
systems. 

Figures IT^ a -IT^d represent the X10Y3 cluster at four 
points in Fig. ^2 with XioY3(cr, e) coordinates (0.8,0.5), 
(0.8,1.0), (0.8,1.5) and (0.8,2.0), respectively. The num- 
ber of inherent structures available to the X10Y3 cluster 
in all four cases is more than 6000. We show only the 
lowest 200 inherent structures. The global minimum of 
each system is labeled by the number 1 and contains 
as a recognizable component the core structure shown 
in Fig. El9. In Fig. Ela inherent structure 1 is linked 
to inherent structure 6 with the isomerization barrier 
which value is AEi^6=0-779exx- The core structure cor- 
responding to inherent structure 6 is different from the 
core structure corresponding to the global minimum. In 
Fig. ^Ib inherent structure 1 is connected to inherent 
structures 7 and 9 whose energy values are -43.177exx 
and -43.034exx, respectively. Both inherent structures 
7 and 9 contain core structures that are different from 
each other and from the one associated with inherent 
structure 1. The isomerization barriers between inherent 
structures 1 and 7 and 1 and 9 are AEi^7=2.138exx and 
AEi^9=2.252exx, respectively. Figures [Hlc and Fig. [Hid 
show that increasing of the value of e increases isomer- 
ization barriers that link inherent structure 1 with inher- 
ent structures 9 and 12, respectively. Numerically, these 
barriers are AEi^9=2.630exx and AEi^i2=2.883exx, re- 
spectively. It can be seen from Fig^^c that by increasing 
the height of isomerization barriers the lifetime can in- 
crease by twenty-one and ten orders of magnitude at 5 
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(12.1) (12.2) 




(12.3) (12.4) 




(12.5) (12.9) 



FIG. 12: Plots of selected X10Y3 structures for various (cr, e) 
values identified in Fig. The number of the structures 
correspond to the regions labeled in Fig. [TTl The core struc- 
ture for the system labeled by (12.9) is not a stable energy 
structure of the bare Xio system. 

and 10 K, respectively. These results are specific demon- 
stration of goal (3) stated earlier. 



IV. CONCLUSIONS 

As stated at the outset, the general theme of the 
present work is to explore the extent to which one can 
induce controllable structural modifications in clusters. 
One reason for pursuing such a development is the possi- 
bility that such modifications might be a general tech- 
nique for producing materials that have "interesting" 
properties (electronic, magnetic, optical, thermal, etc.). 
Depending on the application, one could envision such 
clusters being of use either directly, or, if they could be 
made sufficiently robust, as precursors in subsequent as- 
sembly of yet more complex materials. 

Both current and previous work^^ indicate that such 
controllable modifications are possible in model Lennard- 
Jones systems. In particular, we have demonstrated that 
by introducing impurity atoms of varying size and in- 
teraction energies, we can produce core-atom conform- 
ers that correspond to a variety of non-minimum energy 
homogeneous isomers. We have also shown that it is 
also possible to use such impurities to generate core-atom 
structures that posses no (stable) homogeneous analogs. 
We have now demonstrated such capabilities for both the 
simple X5Y2 and X7Y3 systems and for the larger, more 
complex X12Y1, X11Y2, and X10Y3 binary clusters. 

Finally, in the present work we have investigated the 
issue of the stability of our modified clusters. Based on 
simple transition-state estimates of the rates of isomer- 
ization of the various clusters, we have demonstrated that 
it is possible both to induce and to stabilize a variety of 
structural modifications. In the companion^^ paper we 
investigate how the structures of the underlying poten- 
tial energy surfaces explored in this work are reflected in 
the thermodynamic properties of the systems. 
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